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In most studies of equivalence principle violation by solar 
system bodies, it is assumed that the ratio of gravitational to 
inertial mass for a given body deviates from unity by a param- 
eter A which is proportional to its gravitational self-energy. 
Here we inquire what experimental constraints can be set on 
A for various solar system objects when this assumption is re- 
laxed. Extending an analysis originally due to Nordtvedt, we 
obtain upper limits on linearly independent combinations of 
A for two or more bodies from Kepler's third law, the position 
of Lagrange libration points, and the phenomenon of orbital 
polarization. Combining our results, we extract numerical 
upper bounds on A for the Sun, Moon, Earth and Jupiter, 
using observational data on their orbits as well as those of the 
Trojan asteroids. These are applied as a test case to the the- 
ory of higher-dimensional (Kaluza-Klein) gravity. The results 
are three to six orders of magnitude stronger than previous 
constraints on the theory, confirming earlier suggestions that 
extra dimensions play a negligible role in solar system dynam- 
ics and reinforcing the value of equivalence principle tests as 
a probe of nonstandard gravitational theories. 
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I. INTRODUCTION 

We investigate the consequences of possible violations 
of the equivalence principle (EP) for solar system bodies 
whose ratio of gravitational mass m g to inertial mass rrii 
is given by 



mg/rrii 



1 



(1.1) 



The study of this problem has a long history stretching 
back to Newton |0,pl . In modern times it has been most 
closely associated with Nordtvedt , who investigated 
the possibility that A is proportional to the body's rela- 
tive gravitational self-energy U 



A = r/U 

U = —^r 



(1.2) 



p(x) p{x') 



d 3 x d 3 x' 



I \ p{x) d 6 x , 



where r\ is a universal constant made up of parametrized 
post-Newtonian (PPN) parameters, 77 = 4/3 — 7 — 3 — 
10£/3 - ai + 2a 2 /3 - 2Ci/3 - C2/3 [§. In the standard 
PPN gauge, 7 is related to the amount of space curva- 
ture per unit mass, (3 characterizes the nonlinearity of 



the theory, £ and ctk refer to possible preferred-location 
and preferred-frame effects, and the allow for possible 
violations of momentum conservation (see jlj] for discus- 
sion) . One has r\ = in standard general relativity, where 
7 = j3 = 1. In 4D scalar-tensor theories, by contrast, 
7= (1 + lu)/(2 + lu) and/3 = 1 + w'/[2(2 + w)(3 + 2lu) 2 }, 
where uj — Lu((f>) is the generalized Brans-Dickc parame- 



ter and 



du>/d(f). 



The relative gravitational self-energy U can be cal- 
culated for most objects in the solar system, subject 
to uncertainties in their mass density profiles. Thus, 
for example, U s ~ — 10~ 5 for the Sun, while Jupiter 
has Uj ~ — 10~ 8 ||. More precise estimates are avail- 
able for the Earth and Moon, U E = -4.64 x 10~ 10 and 
U M = — 1.9 x 10~ n respectively [fLo| . With U known, con- 
straints on EP violation (in theories of this kind) take the 
form of upper limits on the PPN parameter ?y. The lat- 
est such bounds (from lunar laser ranging, or LLR) are 
of order |ry| < 1CT 3 ff-[l|. 

There exist, however, theories of gravity in which EP 
violations are not necessarily related to gravitational self- 
energy. In such theories, analysis of solar system data can 
lead only to upper limits on A for the individual bodies 
involved. These limits will in general differ from object 
to object. They will, however, typically be orders of mag- 
nitude stronger than the above-mentioned experimental 
bound on 77 (since they are not diluted by the small fac- 
tor U). In what follows, we obtain numerical constraints 
on A for the Sun, Moon, Earth and Jupiter, and ap- 
ply these as an example to higher-dimensional (Kaluza- 
Klein) gravity, in which EP violations are related to the 
curvature of the extra part of the spacetime manifold. 



II. EQUIVALENCE PRINCIPLE VIOLATIONS IN 
KALUZA-KLEIN GRAVITY 

Theories of gravity in more than four dimensions arc 
older than general relativity 13 1, and most often associ- 
ated with the names of Kaluza Q and Klein |]| (for 
recent reviews, see |l^ , [l7|l ). The extra dimensions have 
traditionally been assumed to be compact, in order to 
explain their nonappearance in low-energy physics. The 
past few years, however, have witnessed an explosion 
of new interest in nort-compactificd theories of higher- 
dimensional gravity |l8|-|2"o|] . In such theories the dimen- 
sionality of spacetime may in principle manifest itself at 



1 



experimentally accessible energies. We examine here the 
theory of Kaluza-Klein gravity JTsJ , and focus on the pro- 
totypical five-dimensional (5D) case, although the exten- 
sion to higher dimensions is straightforward in principle. 

In 5D Kaluza-Klein gravity, objects such as stars and 
planets are modelled by a static, spherically-symmetric 
analog of the 4D Schwarzschild solution known as the 
soliton metric. This may be written (following |2l|j , but 
switching to nonisotropic form, and defining a = 1/a, 
b = (3/a and M* = 2m) 



dS 2 = A a dt 2 - A- a - b dr 2 
- 2 <^ 2 + sin 2 ^ 2Nl 



A b dy 2 , 



(2.1) 

where y is the fifth coordinate, A(r) = 1 — 2M*/r, 
is a parameter related to the mass of the object at the 
center of the geometry, and the constants a, 6 satisfy a 
consistency relation a 2 + ab + b 2 = 1, or 



a = -6/2 ± (1 - 36 2 /4) 



1/2 



(2.2) 



whic h fo llows from Einstein's field equations. Equa- 
tion (2.1) reduces to the 4D Schwarzschild solution on 
hypersurfaces y = const as 6 — > and a — ► +1. 

Solitons differ in interesting ways from 4D black holes, 
as discussed by many authors ]l6| , p7| , pl| -[2^] . Of most rel- 
evance for our purposes is the fact that a soliton's gravita- 
tional mass (as identified from the asymptotic behaviour 
of g o) is given by m g = aM*, while its inertial mass (as 
obtained using the Landau-Lifshitz energy-momentum 
pseudotensor) turns out to be m, = (a + 6/2)Af». Us- 
ing equations (1.1) and (2.2), it follows that 



A = =fV 2 (!- 36 2 /4) 



-1/2 



(2.3) 



It is important to recognize that 6 is not a universal con- 
stant like rj, but depends in principle on local physics f23| | 
and may vary from soliton to soliton. (Birkhoff's theorem 
in its usual form does not hold in higher-dimensional gen- 
eral relativity ^4j.) To constrain the theory, one there- 
fore hopes to apply as many tests as possible to a given 
astrophysical system. 

The classical tests of general relativity have been 
worked out for the soliton metric. Data from long- 
baseline radio interferometry, ranging to the Mars Viking 
lander, and the perihelion precession of Mercury imply 
that |6| < 0.07 in the solar system |25 26|. In what fol- 
lows, therefore, we keep only terms of second order or 
lower in 6. We also drop the negative roots in equa- 
tion ( |2.2j ) [ie, the positive roots in equation (2.3)], in 
order to eliminate the possibility of negative gravita- 
tional and/or inertial mass, and to ensure that the 4D 
Schwarzschild solution is recovered (for y = const) as 
6^0 Ml . This leaves us with 



A = -6/2 



(2.4) 



which may be compared with equation (1.2) for PPN- 
type theories. In Kaluza-Klein gravity, experimental con- 
straints on EP violation translate directly into upper 



bounds on the metric parameter 6, which characterizes 
the departure from flatness of the fifth dimension in the 
vicinity of the central mass. 



III. THE MODIFIED THIRD LAW OF KEPLER 

Let us now proceed to see what experimental con- 
straints can be placed on the parameter A for solar sys- 
tem bodies. Given any inverse-square law central force 
(F = — «/r 2 ), it may be shown that the two-body prob- 
lem admits as solutions elliptical orbits satisfying 



27ra 3 / 2 , 



(3.1) 



where r is the orbital period, a the semi-major axis, and 
fi = mi i JTi^a /(?n.i,i + m.i )2 ) the reduced mass of the sys- 
tem p7| . (Here m^k is the inertial mass of the kth body.) 
We are concerned in this paper with the case in which 
Gm Si im g! 2, where the gravitational masse s m a ± are 



related to the inerti al o nes m^fc by equation (LI). Re- 
arranging equation ( |3.l[ ), one then obtains the following 
modified form of Kepler's third law 



G(mi + m 2 + m 2 Ai + miA 2 ) 



co 2 a 3 



(3.2) 



where to = 2tt/t is the mean orbital angular frequency, 
and rrik refers (here and elsewhere in the remainder of 
this paper) to the gravitational mass of the kth body. 
Equation ( |3.2| ) was first obtained by Nordtvedt using a 
somewhat different argument ||. 

If all bodies have the same value of A, we will not see 
violations of the EP. The common part of Ai and A 2 can 
be absorbed into a rescaled gravitational consta nt. O ne 
way to see this explicitly is to rewrite equation ( |3.2| ) in 
the form 



G(l + Ai)(mi + m 2 ) + Gmi(A 2 - Ai) 



2„3 



lu a 



(3.3) 



There are two modifications to Kepler's third law here: 
a rescaling of the value of G in the first term, and a 
completely new second term, which depends only on the 
difference A 2 — Ai. This latter term is a clear manifes- 
tation of EP violation in the system. 

Nevertheless, we can put experimental constraints on 
the values of Ai and A 2 , even in case s where they are 
equal. To see this, divide equation (3.3) through by Gm\ 
and rearrange to obtain 



Gtoi 



m? 
1 + — 
mi 



= (A 2 -Ai) 



mi 



(3.4) 



To within experimental uncertainty, of course, the left- 
hand side of this equation vanishes for planetary two- 
body systems; this is a statement of Kepler's third law. 
(More accurately, we may say that in standard theory, 



2 



where Ai and A2 are assumed to vanish, a statistical 
best-fit value is chosen for Gm\ in such a way as to force 
the left-hand side as close to zero as possible for all sys- 
tems observed.) Even when (A2 — Ai) = 0, therefore, 
observation imposes an upper bound on the possible size 
of the last term on the right-hand side. This reflects the 
fact that we have excellent experimental data, not only 
on possible violations of the EP, but also on the value of 
Gm\ (as well as u>,a and m 2 /mi) for most solar system 
objects. (The situation is less satisfactory for G and mi 
considered individually |2^] , but these quantities are not 
needed here.) 

We can use this extra knowledge to set limits on the 
sum, rather than the difference of Ai and A2. In combi- 
nation with more purely EP-based constraints from the 
three-body problem, this will allow us to extract limits 
on the individual terms Ai and A2. We therefore collect 
terms in Ai and A2 and recast equation (3.4) in a form 
more suitable for contact with observation 



m,9 . 

— A! + A 2 = 
mi 



m s 
mi 

1 



f- 

m 2 
mi 



(3.5) 



where k is the Gaussian constant, a the semi-major 
axis of the orbit, and A the length of the astronomi- 
cal unit (AU). The quantities k and A are related via 
k 2 A 3 = Gm s to the gravitational mass of the Sun 
(which is actually inferred in practice from a statistical 
fit to experimental data on A). The length of the AU 
(A = 1.5 x 10 11 m) is currently known from planetary 
ranging data to better than ±20 m pg| . 

To within exp erimental accuracy, the right-hand side 
of equation ( p.5|) vanishes (as remarked above; Kepler's 
third law). Our upper limits on the left-hand side are 
then just the sum of the relevant uncertainties in m s /mi, 
to, a, A and m 2 /mi (k is an exact quantity). 

For the Sun- Jupiter system (mi = m s and 1112 = ^1.1), 
uj{= 0.53 rad/yr) has an uncertainty which may be con- 
servatively estimated at 1 arcsec/250 yr Q. The mean 
Sun- Jupiter distance {— 5.2 AU) is likewise uncertain to 
about 200 km, based on post-fit residuals to all observa- 
tional data (not individual ranging measurements to 
Jupiter, which are good to as little as 3 km |3l|]). And 
the mass ratio m s /m,j(= 1047) is known to ±0.0005 [p9| . 
Combining these numbers, we find that 



A s 
1047 



+ A. T | < 8 x 10~ 7 



(3.6) 



The bulk of this is due to the uncertainty in a, with a 
small remainder coming from uncertainty in lu. The un- 
certainties in A (that is, in Gm s ) and mj/m s are unim- 
portant by comparison. 



For the Earth-Moon system (mi 



m E ,m 2 



m. 



uj{= 0.23 rad/day) has an uncertainty of less than 
0.03 arcsec/century [|52|. The mass ratio m s /(m E + 
m M )(= 3.3 x 10 5 ) is known to better than ±0.002 @, 



while m E /m M (= 81) is good to ±0.00002 @. The mean 
Earth-Moon distance (— 3.8 x 10 5 km) has been estab- 
lished to within 0.8 m by LLR data |33|. Let us take the 
uncertainty in the length of the semi-major axis a to be 
less than twice this figure, or 1.6 m. We then obtain 



1 81 



< 1 x 10" 



(3.7) 



In this case, while uncertainty in a is still responsible 
for about three-quarters of the total, there are also sub- 
stantial contributions from uncertainties in m s /m E and 
m M /m E . The decrease in importance of error in a (rel- 
ative to the Sun-Jupiter case) may be attributed to the 
high quality of the LLR data. 



IV. MIGRATION OF THE STABLE LAGRANGE 
POINTS AND THE TROJAN ASTEROIDS 

We turn next to the effects of EP violation on three- 
body motion in the solar system, beginning with the 
problem in which a test body is located at one of the 
stable Lagrange points (L = L± or L5) in the orbit of 
one massive body M 2 (a planet or moon, say) around 
another (M\ ) . Nordtvedt || has shown that the location 
of L is sensitive to Ai ^ for the largest of the three 
masses. We follow his notation, which is based on the 
assumption of circular orbits, but extend the calculation 
to cases in which all three o bjec ts can in principle display 
EP violations of the form (|l.l[) . 

The point L is located at points where radial accel- 
eration toward the M1-M2 center of mass (C, say) is 
precisely balanced by centrifugal acceleration, and where 
there is no net acceleration perpendicular to this direc- 
tion. The situation is illustrated in Fig. 1, where r\ and 
r 2 denote the distances from L to Mi and M 2 respec- 
tively, while r, Ri and R 2 are the distances of L, Mi and 
M2 from C. It is straightforward to show || that 



(ri/r 2 ) 3 -(l + Ai)/(l + A 2 ) 



(4.1) 



The effect of the A-terms, in other words, is to "de- 
equilateralize" the triangle made up by L, Mi and M 2 . 

It is further possible to solve explicitly for n and r 2 ; 
one finds that 



[ri/(Ri+R2)} 3 = (1 + A L )/(1 + A 2 ) 
[r 2 /(Ri+R 2 )f = (l + A L )/(l±Ai) 



(4.2) 



Denoting the unperturbed values of r\ and r 2 by R (= 
i?i + R 2 ), and assuming that the A-terms are much less 
than unity, we see that they lead to small movements of 
L toward (or away from) Mi and M 2 , as follows 

«5ri«-(A 2 -A L )(-Ri+-R 2 )/3 
5r 3 w-(Ai-A L )(iJi+i?2)/3. 



In the limits A L — > and A2 
to those given in || . 



(4.3) 

0, these results reduce 
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Experimentally, one could use these results to obtain 
constraints on (A& — A L ) by ranging to a satellite at L and 
looking for nonzero values of Sr^ (k — 1 or 2). In the case 
of the Earth-Moon system (Mi = Earth, M 2 = Moon) 
this would require an artificial satellite such as the SOHO 
spacecraft (which is however located at L\ rather than L4 
or L5). In the case of the Sun- Jupiter system (Mi = Sun, 
M2 = Jupiter) one could instead make use of the larger 
Trojan asteroids, of which 413 have now been identified 
|Pq| , and some 2 x 10 6 estimated to exist with radii over 
1 km p6[ | . Such a procedure would however be hindered 
by low signal-to-noise levels and uncertainties involving 
factors such as surface topography. 

An observational quantity offering greater promise is 
the angular position 9 of the stable Lagrange point rel- 
ative to the line through Mi and M2 (see Fig. 1). For 
different values of r% and r 2 , one will measure different 
angles 89. From the law of sines 



tan(0 + 59) 



(R-5r x ) wa.fi 



(4.4) 



Ri - (R-5r x ) cos fi ' 
where fi may be determined from the law of cosines 

[R-5r 2 f = R 2 + {R-5r x ) 2 

-2R(R -5r x ) cos fi . (4.5) 

Differentiating, and linearizing in the small parameters 
5ri and 5r 2 , one can show that 



(2gi + R 2 ) Sri - (Ri + 2R 2 ) Sr 2 
V3 (R\ + R X R 2 + R\) 



(4.6) 59 



Inserting equations (4.3) into this expression, we obtain 
the general result 



59 



(Ri + R 2 



[(R X +2R 2 )(A X -A L 



3V^ (Ri + R1R2 + R\ 

-(2R 1 + R 2 )(A 2 - A L )] . (4.7) 



That this is a pure EP-violating effect can be seen from 
the fact that only differences of A-terms appear. In cases 
where Ai = A2 = A L , there is no migration of the La- 
grange point and 59 = 0. 

Let us now specialize to cases in which Mi is much 
more massive than M2, so that R 2 ^> R x . (This is not a 
serious limitation as it describes all situations of interest 
in the solar system.) We also follow Nordtvedt || in 
assuming that the test object at L differs from the much 
larger masses Mi and M 2 in such a way that A L <C A x 
and A L <C A2. This leaves us with 



1 



3-y/3 



(2Ai-A 2 ) 



(4.8) 



a result that can be used to set constraints on (2Ai — A2) 
if we have experimental data on any possible angular shift 
of the Lagrange points from their expected locations 60 
degrees behind and ahead of Jupiter. 



It happens that we do have good observations of a 
number of Trojan asteroids, the first of which was discov- 
ered over ninety years ago. To locate their mean angular 
position with sufficient precision for EP tests is however a 
significant challenge because these objects undergo libra- 
tions about the Lagrange points, with periods that are 
typically a good deal longer — of order 150 yr |57j — than 
the timescale over which they have been observed. Only 
twelve Trojans (588 Achilles, 911 Agamenon, 1404 Ajax, 
1173 Anchises, 1172 Aneas, 1437 Diomedes, 624 Hector, 
659 Nestor, 1143 Odysseus, 617 Patroclus, 884 Priamus 
and 1208 Troilus) have been observed for sixty years or 
more. So one has to fit orbits to an incomplete arc of 
observations. Moreover the older observational data will 
be subject to larger random scatter than more recent 
measurements. There are also several potential sources 
of systematic error, including observational selection ef- 
fects and the nonuniform distribution of asteroids, which 
may not necessarily cancel themselves out as the libration 
centers follow Jupiter around the sky. 

The problem has nevertheless been seriously tackled 
in a pair of papers by Orellana and Vucetich p7| , p8[ . 
(The authors set out to obtain constraints on the PPN 
parameter 77, but we may relate this simply to 59 by 
a constant factor, given as 59 = 1.26^ arcsec in Q.) 
A formal least-squares statistical fit (including a model 
for systematic errors) produces the rather surprising re- 
sult 59 = —0.21 ± 0.09 arcsec, which could perhaps be 
taken as a potential EP-violating signal. A more con- 
servative procedure, however, leads to a final value of 
-0.18 ±0.15 arcsec, which is consistent with no sig- 
nal. Based on this latter limit we adopt the upper bound 
\59\ < 0.33 arcsec. This compares reasonably with plane- 
tary angular measurements, which are typically accurate 
to between 0.1 and 1 arcsec |e3|. In combination with 



equation (4.8) this bound leads to the constraint 



(4.9) 



|2A S - A., I < 8 x 10" 



This is an order of magnitude weaker than the limit ( |3.6| ) 
on the sum A s /1047+A., for the Sun- Jupiter system from 
Kepler's third law, and almost three o rde rs of magnitude 
weaker than the corresponding limit (3.7) for the Earth- 
Moon system. Nevertheless equation (i.£) will prove use- 
ful in constraining the value of Aj. 



V. ORBITAL POLARIZATION AND THE LUNAR 
NORDTVEDT EFFECT 

We turn next to the best-known probe of EP viola- 
tions in the solar system, the polarization of a two-body 
orbit in the field of a third body, also known variously as 
the Nordtvedt Effect and the Lunar Eotvos experiment. 
This was investigated by Nordtvedt in the context of the 
Earth-Moon-Sun system and the Sun-Earth-Jupiter 
system \ t], and has been the subject of numerous studies 
since ggfj^y]. We follow the approach of Will @, but 
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extend his notation to general three-body systems. The 
situation is illustrated in Fig. 2. 

Including their mutual attraction (but neglecting any 
quadrupole or higher moments), the tw o massive bodies 
Mi and M 2 described by equation (O) will fall toward 
M3 with different accelerations given by 



ai = -(1 + Ai)(m 3 xi/rf - m 2 x/r 3 ) , 
a 2 = -(1 + A 2 )(m 3 x 2 /r^ + m lX /r 3 ) , 



(5.1) 



where Xi and x 2 are the 3-vectors connecting M 3 to Mi 
and M 2 respectively, and x is the 3-vector from Mi to 
M 2 . (The lengths of these vectors are 7*1, r 2 and r respec- 
tively.) The relative acceleration a = a 2 — ai between Mi 
and M 2 may then be written 



a = — 



X Xi 

//i*-? + (Ai - A 2 ) 7713-3- 



-(1 + A 2 )m 3 [^-? 



(5.2) 



where to* = mi(l + A 2 ) + m 2 (l + Ai). The first term 
is the standard Newtonian acceleration, which gives rise 
to the circular (unperturbed) orbit of M 2 around Mi, 
assumed here to be coplanar with M3. The second, 
"Nordtvedt term" will be treated as a perturbation of 
this orbit. (We follow Will § in neglecting the third, 
tidal term at this stage.) The acceleration a = d 2 x/dt 2 
and angular momentum h = x x (dx./dt) satisfy 



d 2 r/dt 2 = x ■ a/r + h 2 /r 3 , 
dh/dt = x x a . 



(5.3) 



Defining 8a = (Ai — A 2 )m 3 /ri, linearizing about the 
unperturbed orbit (so that r = 7*0 + Sr and h = ho + Sh) 
and solving the resulting pair of differential equations, we 
find the following periodic oscillation in distance between 
Mi and M 2 

Sr = (Ai - A 2 ) A EP cosD , (5.4) 
where A Er , the EP-violating amplitude, is given by 



A E p = 



1 + 2w 2 /(cj 2 — uj\) 
2(w 2 /o;i) -1 



n 



(5.5) 



and D = (u 2 — u)\)t is the synodic phase. Here u>i is 
the orbital angular frequency of Mi about M 3 , w 2 is the 
orbital angular frequency of M 2 about Mi , and we have 
made use of the relations m*/rQ 



m 3 /r 2 = wfri 



and ho/r 2 , = L02 



Equation ( p.4|) represents a polarization, or alignment 
of the orbit of M 2 about Mi along the direction either 
toward (if Ai > A 2 ) or away from M3 (if Ai < A 2 ). 
For the Earth-Moon system (Mi = Earth, M 2 = Moon, 
M 3 = Sun), we use lu± = 1.991 x 10~ 7 rad/s, w 2 = 
2.662 x 1(T 6 rad/s and n = 1.496 x 10 11 m to obtain 
A EP = 1.84 x 10 10 m. Due to amplification by tidal ef- 
fects, however, this figure is thought to be a consider- 
able underestimate More refined calculations, both 



analytical and numerical, now indicate a final value of 
A EF = 2.9 x 10 10 m UfiHOJ. 

All these amplitudes are, of course, far larger than 
any observed anomalous fluctuations in the Earth-Moon 
distance at this frequency. (The largest known term at 
the synodic frequency has an amplitude of approximately 
100 km, but can be modelled so accurately that it does 
not contribute to the uncertainty in Sr.) Indeed, anal- 
ysis of LLR data now constrains any such fluctuations 
to be less than 1.3 cm in size pHll||. Using the latest 
above-mentioned value for A EP , we infer an upper bound 



I A K — A M I <4.4x 10 



-13 



(5.6) 



For PPN-type theories, in which A E = i]U E and A M = 
r]U M with the values of U E and U M given i n S [| this 
leads to the constraint \tj\ < 1.0 x 10~ 3 jljlij, which 
is currently the strongest limit on this parameter. (It 
has been weakened only slightly, to \rj\ < 1.3 x 10~ 3 , 
by a recent study of possible masking by composition- 
dependent effects fl2|.) 



VI. 



ORBITAL POLARIZATION AND THE 
SOLAR NORDTVEDT EFFECT 



One can also apply equation (5.4) to the Sun-Earth 
system with Jupiter as the third mass (Mi = Sun, 
M 2 = Earth, M 3 = Jupiter), as suggested in ||]. This 
might be dubbed a solar Eotvos experiment since it is 
the difference in accelerations of the Earth and Sun, and 
the size of anomalous variations in the solar-terrestrial 
distance, that are of interest. For this case we have 
= 1.679 x 10~ 8 rad/s, u 2 = 1.991 x 10~ 7 rad/s 
m, 



and r\ — 7.783 x 10 11 m, resulting in an uncorrected 



Nordtvedt factor of A EP = 1.09 x 10 11 m. 

A realistic upper limit on anomalous fluctuations in 
the Sun-Earth distance at the relevant frequency (corre- 
sponding to a period of 1.09 yr) might be several times 
the uncertainty in the length of the AU (±20 m); we take 
here a value of 100 m. (There are ordinary classical per- 
turbations of amplitude ~ 2000 km at this frequency || , 
but these can be compensated for very accurately, since 
the ratio of Jovian to solar mass is known to better than 
a part in 10 6 p9|.) Using equations (5.4) and (5.5), we 
then obtain 



I A= - A R I < 9 x 10" 



(6.1) 



which is some four orders of magnitude stronger than 
the constraint on A s and A,, obtained from the Trojan 
asteroids in equation ( f4.9| ). This is a reflection of two 
things: the uncertainties discussed in § [V in connection 
with the motion of these bodies, and the fact that our 
measurements of the distance to the Sun, like that to the 
Moon, rest ultimately on the great precision of ranging 
data (in this case, from the Viking lander rather than the 
Apollo retroflectors) . 
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To complete the set of constraints, we mus t fin ally ap- 
ply the modified Kepler's third law, equation (|3.5| ), to our 
last pair of bodies. For the Sun-Earth system (mi = m s 
and mi — m E ), uncertainty in u>(= 2t: rad/yr) will be 
insignificant. For definitcness one could perhaps use the 
total drift error in the Earth's inertial mean longitude, 
which is 0.003 arcsec/century psfl . As in the Earth- 
Moon case, let us take the uncertainty in semi-major 
axis a (here 1 AU) to be 40 m, or twice that in the 
length of the AU itself. Uncertainty in the mass ratio 
m E /m s (— 1/334,000) may be obtained from the limits 
£[m s /(m E + m M )) < 0.002 (33) and 5(m E /m M ) < 0.00002 
p4j , giving a figure of S(m E /m s ) = 2xl0~ 14 . Combining 
these numbers, we find that 



A. 



1 334, 000 



AJ < 7 x 10~ 9 



(6.2) 



This is due almost entirely to the uncertainty in m E /m s 
— small though that is — with the tiny remainder com- 
ing from uncertainty in the values of a and A. 



VII. INDEPENDENT LIMITS ON A 

We now proceed to extract individual const rain ts on 
the parameters A& by combining equations (3.6) and 
(p~9)), ( |3~7| ) and <^), and (|]|) and (f~2]) respectively. 
All three pairs of inequalities may be written in the form 



|Ai +ciA 2 | < e x , 

|Ai -C3A2I < £2 , 



(7.1) 



where ci,C2,ei and £2 are positive. We wish to obtain 
upper limits on the qu anti ties A^. 

Squaring equations (7J), taking their sum and differ- 
ence, and solving the resulting pair of quadratic equa- 
tions, we find that 



|Ar|< 



|A 2 | < 



|(ci + c 2 )(ei T £2) - (ci - c 2 )(ei ± e 2 )| 



2(ci + c 2 ) 



|ei zbe 2 | 
ci + c 2 



(7.2) 



where the signs must be evaluated together (ie, if the 
upper sign is selected for Ai, then it must be selected for 
A 2 as well). 

We apply these results to the Sun-Jupiter system 
(l=S,2=j) by substituting ci = 1047 and ei = 8 x 10~ 4 
from equation ( |3.6| ) (Kepler's modified thi rd la w), and 
C2 = 0.5 and e 2 = 4 x 10~ 6 from equation (L9) (migra- 
tion of the stable Lagrange point) to obtain 



|A S | < 5 x HT 6 , 
I A., I < 8 x 10" 7 . 



(7.3) 
(7.4) 



Simil arly, for the Earth-Moon system (1=e,2=m), equa- 
tion C7|) gives ci = 81 and e x = 8 x 10~ 7 (Kepler's 



modified third law), while equation ( J5.6| ) gives c 2 = 1 
and e 2 = 4 x 10~ 13 (lunar Nordtvedt effect), so that 



|A E | < 1 x 10" 
IAJ < 1 x 10" 



(7.5) 
(7.6) 



Finally, for the Sun-Earth system (1=S,2=e), we take 
ci = 334,000 and e x = 2 x 10~ 3 from equation (|6~J) 
(Kepler's law), while c 2 = 1 and 62 = 9 x 10~ 10 from 
equation ([O]) (solar Nordtvedt effect), so that 



|A S | < 8 x 10~ 9 
|A E | < 7 x 10~ 9 



(7.7) 
(7.8) 



The limit ( |7.7|) obtained for the Sun from the solar 
Nordtvedt effect is nearly three orders of magnitude 
stronger t han that derived from the Trojan asteroids, 
equation ( |7.3| ). This is perhaps not surprising, given 
the discussion in § IV above. But the fact that the lu- 



nar and solar Nordtvedt effects lead to equally strong 
cons traints on A E for the Earth — equations (7.5) and 

(7.8) respectively — is somewhat unexpected, in view of 
the widespread belief that LLR data provides by far the 
strongest probe of EP violations in the solar system. This 
may be partly explained by the fact that the strength 
of both Nordtvedt effect-based limits ultimately derives 
from ranging measurements, as discussed in § |VT| . 

Statistical effects may also play a role here. As uncer- 
tainties in our orbital and other parameters (x, say), we 
have in each case used residuals from published fits to a 
fixed number of solution parameters. These fits, however, 
do not generally incorporate all the % (they are typically 
sensitive to at most a single EP- violating parameter, rj). 
We have, in other words, relied on more degrees of free- 
dom than are actually present in the solutions. This is 
not necessarily a problem, but will tend to underestimate 
our uncertainties. The results least affected will b e th ose 
based on the lunar Nordtvedt effect — equation ( |5.6| ) - 
for which EP-violating terms have been incl uded i n th e 
solution sets. Our other results — equations (|3.6| ), (0), 

(4.9) , (3.1) and (3.2) — maybe less robust in comparison. 
A fully rigorous future treatment would rely on new sta- 
tistical fits incorporating an independent parameter Afc 
for each body, rather than parametrizing EP violations 
with a single universal constant. 



APPLICATION TO KALUZA-KLEIN 
GRAVITY 



VIII. 



Let us now see what equations ( |7.3| ) ~ ( |7.8[ ) imply for 
5D Kaluza-Klein gravity, in which [as we recall from 
equation ( |2.4| )] A w —6/2, where b is a free parameter 
which characterizes the departure from flatness of the 
fifth dimension in the vicinity of the massive body. Our 
constraints on this parameter follow immediately from 
equations (f7~J), (\T(\), fl7~7| ) and (|7~~ 
marized together as follows 



and may be sum- 
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|6 S | < 2 x 10~ 8 
|6j| < 2 x 1(T 6 
|6 E | < 2 x 1(T 8 
\b M \ < 2 x 1CT 8 



(8.1) 
(8.2) 
(8.3) 
(8.4) 



These upper bounds on b are far more stringent than any 
obtai ned so far by other means. In particular, the re- 
sult (3.3) for the Earth is more than three orders of mag- 
nitude stronger than that which may be expected using 
data on geodetic precession from the up com ing Gravity 
Probe B satellite |^|{||9). Equation (fj) for Jupiter 
is more stringent by five orders of magnitude than the 
limits which may be obtained using light defle ction by 
the giant planet And the upper limit (8.1) for 

the Sun is more than six orders of magnitude tighter 
than those set so far from the classical tests of general 
relativity (light deflection using long-baseline radio in- 
terferometry, Mercury's perihelion precession, and time 
delay using the Mars Viking lander p5|,p6[). 



IX. CONCLUSIONS 

We have looked for the constraints imposed by solar 
system data on theories in which the ratio of gravita- 
tional to inertial mass differs from unity by some factor 
A which may in principle differ from body to body For 
two objects characterized by Ai and A2, upper bounds 
on the sum |Ai + C1A2I have been found from Kepler's 
third law, while the difference |Ai — C2A2I can be con- 
strained by data on the position of Lagrange libration 
points, and orbital polarization in the field of a third 
body (the Nordtvedt effect). (Here c\ and C2 are known 
constants.) Combining these results, we have extracted 
independent upper limits on A for the Sun, Moon, Earth 
and Jupiter, using experimental data on their orbits as 
well as those of the Trojan asteroids. In particular, we 
find that A < 10~ 8 for the Sun, Earth and Moon, and 
1CP 6 for Jupiter. 

As a test case, we have applied these results to five- 
dimensional Kaluza-Klein gravity, in which A can be 
shown to depend on a parameter b of the five-dimensional 
metric which characterizes the curvature of the extra di- 
mension near the central mass. Our upper bounds on this 
parameter are three to six orders of magnitude stronger 
than existing limits on b from the classical (and geodetic 
precession) tests of general relativity, and confirm earlier 
conclusions |l(||26| that a fifth dimension, if any, plays no 
significant role in the dynamics of the solar system. 
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FIGURE CAPTIONS 

FIG. 1. Effect of EP violations on the stable three- 
body configuration of Lagrange. 



FIG. 2. The Nordtvedt effect: two bodies M x and M 2 
fall toward M 3 with different accelerations. 
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